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Abstract 



We propose a single particle boson mechanism of High Tc Superconductivity 
(HTCS) and pseudogap regime. Bosons appear in it due to the coupling of spins 
of the two-dimensional (2D) fermions with statistical magnetic field induced by 
anyon vector potential. The ground state of 2D gas is pure bosonic if gas is 
not dense. At the dense limit of gas the interaction of effective (coupled with 
the statistical magnetic field) spins of bosons leads to the increasing of their 
fluctuations, which destroy the coupling. An experimental phase diagram of the 
hole doped superconducting cuprates discussed in the paper of Tallon and Lo- 
ram might qualitatively and quantitatively be clarified in the framework of this 
mechanism. The vicinity of the structural phase transition to superconducting 
state might strengthen the possible quadratic striction in the sample and the 
phase transition of bosons into Bose- Einstein condensate (BEG), which is re- 
sponsible for the superconductivity (SC), is not second order, but first, close 
to second one. According this treatment the pseudogap regime is the region of 
meta stable bosons, which are out of the BEG. At the pseudogap boundary. Eg, 
the bosons finally undergo the phase transition into fermions. Non-Fermi liquid 
like property of quasi-particles discussed in the literature might be related to 
bosons with spins in the pseudogap regime. 

1 Introduction 

A mechanism of SC based on the BEG of bosons, irrespective of nature of parti- 
cles, has been proposed by Ogg, Ref. [l],in 1946 and then by Schafroth together 
with collaborators, Refs. [21 [3], (see also the review of Ginzburg, Ref. [3| ) during 
the decade since 1950. It seems a single particle boson nature of SG is mysteri- 
ous, because no reason for these particles to appear in the solid materials. Only 
astrophysical objects, like stars, allow [4] at the extremal physical conditions this 
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scenario of SC. It is believed [5] that, as metals at low temperatures, the HTCS 
materials are described by Cooper like pairs of fermions and the set of 2D supercon- 
ducting planes of atoms Cu02 plays the important role. For these 2D systems has 
been suggested the alternative semion (anyons with fractional parameter u = 1/2) 
picture of HTCS, which was latter not supported in the experiment ( the references 
on the theoretical and experimental papers with respect to the subject see in Ref. 
[6]). The extended discussion of this treatment, where anyons are the quasi-particle 
excitations, is outlined in |7j. 

In general, the concept of anyons, as result of non-relativistic Chern-Simons 
Quantum Field theory |6| , discovers the physical richness of two space dimensions. 
The topology of configuration space for orbital motion of particles in 2D allows 
for fractional exchange statistics [8], characterized by a continuous parameter v 
that may attain values between (for bosons) and 1 (for fermions). Particles with 
< < 1 are generically called anyons p]. 

Assuming anyons are spinless, it is believed [6] that particles in this concept can 
continuously transform from canonical bosons into canonical fermions [10]. This 
means, for instance, the symmetric wave function of boson ground state can contin- 
uously transform into antisymmetric one of fermion ground state. In the paper we 
follow this notion. 

The interesting and not yet considered in the literature on anyons problem is the 
relation of fractional statistics and real spins of particles. From standard courses 
of Quantum Electrodynamics (see, for example, Ref. [11]) it is well known that 
particles with the integer number of h/2 spins possess a Fermi statistics with Pauli 
exclusion principle for occupation of one quantum state by particle and antisymme- 
try constraint for the many-body wave function. In the same time, for 2D systems 
the concept of anyons (at v = 1) gives an opportunity to introduce this antisymme- 
try property of the wave function into Hamiltonian. In the paper we assume that 
anyons (like electrons or holes) have the spin h/2. 

The goal of the present paper is to outline the results of paper [12], where we 
have studied the simultaneous effect of spins and fractional statistics and found that 
the value h/2 is crucial for 2D systems. We have introduced in the Hamiltonian of 
anyon gas the Zeeman term of the interaction of spins of particles with magnetic 
field induced by anyon vector potential, i.e. the statistical magnetic field [6l|T3], and 
showed that the calculation of an expectation value for ground state energy exhibits 
the total cancellation of terms connected with statistics. As the cancellation occurs 
at any v this would mean the bosonization of anyons and, at particular case 
= 1, of 2D fermions in the ground state. Expecting that this effect would be 
general for any 2D gas of fermions, we have applied it for the interpretation of 
phase diagram of HTSC and pseudogap regime suggested by Tallon and Loram in 
the review of experimental papers, Ref. |14j . 

Previously, we have derived an approximate analytic expression for the ground 
state energy of charged anyons confined in a 2D harmonic potential [15]. This was 
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achieved by using the bosonic representation of anyons and a gauge vector potential 
to account for the fractional statistics, which allowed working with the product 
ansatz for the N -body wave function. A variational principle has been applied 
by constructing this wave function from single-particle gaussians of variable shape. 
As in many other perturbative treatments of anyons in an oscillator potential (see 
references in Ref . |15] ) our expression for the ground state energy had a logarithmic 
divergence connected with a cut-off parameter for the interparticle distance. Making 
use of the physical argument (see Ref. [6j) that for 7^ this distance has to have 
some finite value, we have regularized the formula obtained for the ground state 
energy by an appropriate procedure that takes into account the numerical results 
for electrons in quantum dots in the case with Coulomb interaction. 

In our treatment, Ref. [12], the normal and superconducting states were sepa- 
rated by gap - the difference of the ground state energies of fermions and bosons 
and thus, at the evaluation of phase diagram for BEC we have used the approximate 
analytic expressions for these energies obtained in our recent paper [16] for the 2D 
homogeneous Coulomb Fermi and Bose gases. 

The paper is organized as follows. After brief discussion in the Section 2 of 
our approach and results [15] for the confined in harmonic potential anyons with 
Coulomb interaction, in the Section 3 we describe the harmonic potential regulariza- 
tion procedure [16j to get the thermodynamic limit and then obtain the approximate 
analytic expression for the ground state energy of the 2D homogeneous Coulomb 
anyon gas. The results of paper [12] about the implicit anyon or sigle particle boson 
mechanism of HTCS and pseudogap regime will be outlined in the Section 4 and we 
summarize and conclude the present paper by the Section 5. 

2 The Coulomb interacting anyons in a 2D harmonic 
potential 

The Hamiltonian of N spinless anyons of mass M and charge e confined to a 2D 
harmonic potential, interacting through Coulomb repulsions, is given by 



Here and pk represent the position and momentum operators of the kth. anyon in 
two space dimensions. 



normal to the 2D plane. The factor u determines the fractional statistics of the 
anyon: it varies between u = (bosons) and u = 1 (fermions). 




(1) 
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We employ a variational scheme by minimizing the expression for the total energy 

^ ^ J^UR)H^t{R) dR 
J^*^{R)^t{R) dR 

with a trial wave function ^xiR) depending on the configuration R = {fi....rjv} of 
the N anyons. 

It is reasonable, in the bosonic representation of anyons when the many-body 
wave function takes the product form 

AT 

^TiR) = n V-TCrl) (4) 
k=l 

to adopt the single-particle trial functions ipTi^k) in the form 

V^T(ffc) = Cexp(^-(a' + i.)M±^^ . (5) 

Here C is a normalization constant and a' a variational parameter. We identify L 
with the characteristic length {fi/MujoY^"^ of the harmonic oscillator. 

When energies are expressed in units of Tlloq and lengths in units of L the nor- 
malized trial wave function reads 



fc=l 



where a = a' + v. 

In evaluating the expectation value E (Eq. ([3])) it is convenient to consider the 
local energy El{R) = ^^\R)H'^t{R) M- In general El{R) is a complex function 

El{R) = RcEl (R) + ilmEL (R) (7) 

with 

N 

ImELiR) = -aJ2iiMrk) + ei'exi(rfc)/c) • n) ■ (8) 



fc=i 



However, evaluation of the expectation value E = f ^j^iR) EL{R)^TiR) dR imme- 
diately yields 

J ^t{R) ImEL{R)-^T{R) dR = , (9) 

and, therefore, the only quantity to consider in the following is ReEL{R). Before 
proceeding, we would like to emphasize that the absolute ground state of the anyon 
system is a non-analytic function of u. Our calculations will simply provide a smooth 
interpolation. 
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In the non-interacting case the local energy is 

^ ^2 I „,2 ,,2 



ReEUR) = E[« + ^^(1 - a') + ^{Mn)?] ■ (10) 



k=l 



The calculation [15j of the expectation value for ReE^^R) gives 



E = ^(m a + -] (11) 



2 V a 
with 

M=l + v\N - l)[ln - G{N - 2)] , (12) 

dE 

and G = 3^/^ ln(4/3), which attains a minimum (— — = 0) for 

da 

ao = N-^/^ . (13) 
Thus, the resulting expression for the ground state energy is 

Eq = N _ (14) 

The logarithmic divergence displayed in £"0 when 5 — > has also been found in 
other approximate perturbative treatments of the problem and is widely discussed 
in the literature (see, for example, Refs. [20l [211 [22] )• Here we have assumed as 
in Ref. [6j that the cut-off parameter 5 cannot be zero for > 0, away from the 
bosonic limit, since it corresponds to the square of the nearest distance between 
the particles (in Ref. yiSj we have supposed that rg ~ L, where tq is the mean 
distance between particles). Thus, for anyons in the parabolic confining potential 5 
is definitely smaller than 1 (in units of L^). 

Wu [17] has computed the ground state energy of N anyons in a 2D harmonic 
potential near the bosonic limit z/ ~ and obtained 

E^ [N + N[N -l)u/2\ . (15) 

To regularize the expression for Eq we have made use of this result by expanding 
Eq, Eq. (jl4p . for — > and identify the leading term in i^^ with the term linear in 
of Eq. ()15p . with the result 



. 1 
= - exp 



1 + z/G(iV - 2) 



(16) 



With this value of the cut-off parameter the final analytic expression for the ground 
state energy is 

Eq = N[l + u{N - l)]^/^ . (17) 
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This formula for large N is consistent (up to a numerical factor) with the ap- 
proximate expression E ~ v^l'^N'^l'^ of Chitra and Sen [22j calculated perturbatively 
from the bosonic end for v > \ /N . 

It should be noted, that due to this regularization procedure the ground state 
energy obtained, Eq. (fT7|) . is not an upper bound to the ground state as one would 
expect from a variational principle. This is a consequence of the fitting of 6 to the 
result of Wu Ref. [IT], which for = 2 leads to a square root dependence in i/, 
while the exact result for this case gives a linear dependence. On the other hand, 
Eq. (fTTl) applies for the whole range of parameters of the system N, u, and ujq. 

We now include the effect of the Coulomb repulsions between anyons 



L 
2aB 



N 

E 



1 



(18) 



in the expression for the real part of local energy ReEL{R), Eq. (jlOp . 

The Coulomb interaction part contributes with A^(A^ — 1) integrals of the form 

/ ^'r(i?) j-^ — r^riR) dR. These integrals have been evaluated in Ref. [T5| and the 
averaged (real part of the) local energy is 



E 



N 



M a + - + 2M a^/A , 
a J 



with 



M 



2, 



2 an ' 



(19) 



(20) 



dE 



where is the Bohr radius. The extremum condition — — = leads to the equation 

da 



- MX -M=0 

for X = 1/a^/^. The minimum energy is given by the expression 



2 



M ^2 

— 2 + Aq + —— 



and it is achieved at the point 

Xo = {A + Bfl^ + [-{A + B) + 2{A^ - AB + b2)1/2]1/2 

where 



A 
B 



Myi28 + ((AA/12)3 + (>[Vl28)2)^/2' 
"XV128 - ((A/'/12)3 + (7WVl28)2)^/2 



1/3 
1/3 



(21) 



(22) 



(23) 



(24) 
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Figure 1: Coulomb interaction parameter L/as {tq k, L, [15]) dependence of the 
ground state energy for 7-10 electrons calculated by variational [23] and fixed-node 
quantum Monte Carlo methods [24] - dashed curves (results of both calculations are 
indistinguishable in these curves) and by formula (j22p - solid curves. 

Again, the ground state energy Eq, Eq. (j22l) . has a logarithmic divergence in the 
limit 8 ^ {) and the quantity TV should be regularized. 

In order to determine the cut-off parameter 5, and due to the lack of analytic 
results, we had to fit to known numerical results for the ground state energy at 
special values of the parameter L/qb- 

In Figure 1 we have compared the ground state energies calculated for 7-10 
electrons using Eq. (j22p . with the non-interacting J\f = \ -\- v{N — 1), to variational 
[23] and fixed-node quantum Monte Carlo calculations [24J. 

3 Harmonic potential regularization 

In the Section 2 we have outlined a variational procedure for the ground state energy 
of confined interacting anyons (starting from the bosonic end) and achieved, after 
regularization of a logarithmic expression by means of a cut-off parameter for the 
particle-particle interaction, approximate analytic formulas in terms oi N , u and 
L/aB- For the non-interacting anyon system we found 



(25) 
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while for the interacting anyon system the approximate analytic ground state energy 
is given by 



M ^2 2>I 

-^0 ^0 



(26) 



with the expression for Xq, Eq. ([23|) . In these expressions we have used M = 
1 + i/(iV - 1) and Eq. ([20]) for X. 

In order to obtain the corresponding expressions for the 2D homogeneous anyon 
gas, which have been derived in the paper and will be briefly outlined in this 
section, we had to flatten out the parabolic confining potential while increasing the 
number N of anyons, keeping the density p = N/S = 1/vrrg constant, i.e. we per- 
formed the thermodynamic limit while making the confining potential disappearing. 
Here vrrp is the area of the jellium disc carrying the positive countercharge and the 
mean particle distance rg = asrs can be expressed in units of the Bohr radius by 
the dimensionless density parameter r<j. 

Without Coulomb interaction and in the case of fermions {v = 1) the ground 
state energy of the 2D homogeneous electron system of density p is determined by 
the Pauli exclusion principle and is given by 

Eo{p) = Tif^pN/M , (27) 

while from (j25p we have 

EQ{N,v = l) = hujQN'^/^ . (28) 

In the thermodynamic limit both expressions have to become identical and we obtain 
the relation 

cuo(iV) = Trhp/{MN^/^) , (29) 

which means that, in fact, the thermodynamic limit (A^ oo) is obtained for 
vanishing parabolic confinement potential. We can extend this consideration for the 
fermionic limit, u = 1, to the general anyon case, u ^ l,hy assuming instead of (|27p 
the relation 

Eo{p, u) = Trn'^pN(j){u)/M , (30) 

where the function (pii^) is still to be determined under the constraint 0(z^ = !) = !. 
This form is in the accordance with the fact that close to the bosonic limit (v ~ 0) 
the ground state energy of the infinite anyon gas depends linearly on i/ [25\ [26l [27] . 
Consequently, we have for this case 

uo{N, u) = 7Thpf{u)/{MN^/^) (31) 

with another unknown function /(z^) and the constraint /(z^ = 1) = 1. It turns out 
(/)(i^) is determined by /(i^). 

In the thermodynamic limit and including the Coulomb interaction, the 
parabolic confinement has to be replaced by the jellium contribution, which for 
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a disc of radius Rq (containing N counter charges) gives a potential energy contri- 
bution [28] 

Vin) = -pf . (32) 

Js \rk - r\ 

Here S = ttRq is the area of the jehium disc for charges and we have Ro = N^^'^rQ. 
For 1^ = 1, there is a relation between characteristic length L of the oscillator with 
the mean particle distance L = N^^'^r^ and we have found rg ^ L <C i?o for ^ 1. 

In the general case of z/ 7^ 1 and the Coulomb interaction, the approximate 
analytic expression for the ground state (|3Up can be written in the form 

Eo{u, r,) = 7Th^pN4>{u, rs)/M , (33) 

thus, in a similar way to generalize /{v) by f{i^,rs). 

Before deriving the approximate analytic expression for the ground state energy 
of Coulomb interacting anyon gas, in Ref. pS] we have obtained the formula for one 
of the 2D Coulomb Bose gas at high densities. Performing the calculation for 2D 
case, which is the analog of Foldy's one, Ref. |29j, for the 3D case, we have found 

I = -CBcr^'/' , (34) 

where cbg = 1.29355. Thus, we had an exact analytic expression for the ground 
state energy per particle for the 2D Coulomb Bose gas valid at high densities. It 
has next used together with the known expression for the 2D Coulomb gas in the 
low density limit (the 2D Wigner crystal) to derive a form for the unknown function 
/(z/, r<j) introduced above in this section. 

In the paper Ref. [16] we have also obtained the spectrum of collective excita- 
tions of the 2D Coulomb Bose gas. It has the form 

Ek = [TiH^e^pk/M + (n2fc2/(2M))2]i/2 ^ (35) 

which for small k is the 2D plasmon dispersion and approaches for large k the free 
particle dispersion. 

For the derivation of an approximate analytic expression for the ground state 
energy of the Coulomb anyon gas, in paper [16j we have at beginning calculated the 
contribution of the jellium term, Eq. ()32p . into expectation value for energy. We 
have obtained 

pE / r *T(r,) ^^$^(4) = ^iVVV2aV2 . (36) 

y._^JS J -co \Tk ~ f\ ^ 

It was achieved at the condition ^1. A validity of this condition had been 
shown in [16] for almost all numerical values of v and r,. 
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To bring the contribution of the jelhum term in our approximation together with 
the Coulomb interaction evaluated in [l5j and outlined in the previous section we 
might write the ground state energy as in Eq. (j26p but with hujQ = h'^/{ML'^) and 



M 



vr 



OB 



1 



23/2 



1 



(37) 



However, this expression for A4 we have replaced (the reason for that see |16j ) by 
A4 = —cwcN^^'^L/aB, where the constant cwc has been fixed by the ground state 
energy of the Wigner crystal (limit 1). 

For sufficiently large A^, we used M = —cwc^'^^^^s and M = uN to write the 
expression for the ground state energy, Eq. ([26]), with Xq = N^/'^Kx, where 



Kx = (Ka + Kef'^ + [-{Ka + Kb) + 2{Ki - KaKb + 



2Nl/2nl/2 



and 

Ka = [^2/128 + {{u/l2f + (i^Vl28)2)^/^' 
Kb = [i^Vl28 - {{u/l2f + {K"^ /l2Sff'^ 

with K = cwcfs- It took the form ( in Ry units) 



1/3 
1/3 



(38) 



(39) 



2K\ 



+ 



K 
K^ 



(40) 



If we remember that we have introduced a function f{u,rs) in L, therefore, now 
K = cwcTs/ f^^'^{i^^^s) and the final expression for the ground state energy per 
particle is 



2f{v,rs 



2K\ 



+ 



K 



X 



K 
K^ 



(41) 



We have determined the function f{v,rs) by fitting Eq. (j4ip to the known asymp- 
totic limits for the ground state energy of spin polarized electrons and Coulomb 
2D Bose gases at very small and very big values of r^. At very big the ground 
state energy does not depend on statistics and equals the energy of the classical 2D 
Wigner crystal [30], Ey/c = — 2.2122/rs. For the Bose gas at small r<j we used Eq. 

For the bosons at = we had from Eq. (j4ip 



'?o(0,r^ 



2/3 



(42) 



with 



(43) 
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where c^}, = 2.2122 or cwc = 3.2903. 

The asymptotics of the ground state energy, Eq. (j^T|) , for i/ / and at big 
has the form 

The function /(z/, r^) has to fulfil the constraints: f{y = 1, = 0) = 1 for the dense 
ideal Fermi gas; /(i^, = 0) = v^l'^ for the ideal anyon gas close to the bosonic limit 
(see [25l[26l[2Z] ), and /(0,r^) given by Eq. (gSD for the 2D Coulomb Bose gas. The 
form 

/(z., r,) « - r,)e-- + , (45) 

is consistent with these requirements. 

Using this form of f{u,rs) in Eq. (jH]), we have calculated the ground state 
energy per particle in the boson (u = 0) and fermion (zv = 1) limits and compared 
the results for the latter case with the data obtained by Tanatar and Ceperley [31] 
for the spin polarized electron system. 

If we look at the results for large in Figure 2 with ground state energies 
calculated from Eqs. ([^T]) and (jl^ (the results coincide on the plotted scale, solid 
line) and from Eq. (I42p (dash-dotted line) in comparison with the data of [31], we 
see the results from our approximate analytic formula are very close to the "exact" 
data of [31] for the spin-polarized 2D electron system obtained numerically. 



4 Implicit anyon or single particle boson mechanism of 
HTCS and pseudogap regime 

We outline in this section results of paper [12j, where it has been shown that inside 
of the ground state the interaction of spins of anyons with statistical magnetic field 
can induce the bosonization of these particles, i.e. the transformation of anyons (as 
also 2D fermions ) into bosons. 

In the paper [12] we have considered the gas of spinless anyons, which was 
described by the Hamiltonian, Eq. ([1]), with positive background jellium term, Eq. 
(I32p . included inside of the first sum. 

We have introduced in the Hamiltonian the term 

n ^ 

MT.^s.h (46) 
fc=i 

with statistical magnetic field [BJ [13] 



(47) 
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Figure 2: Ground state energy of 2D spin-polarized electrons vs. density parameter 
rs calculated from Eqs. (j4ip or (j44p for u = \ (solid line) in comparison with the 
data of Ref. [31j (dotted line). The dash-dotted line is obtained from Eq. (I42p for 
the ground state energy of the 2D Coulomb Bose gas. 

which can be derived if calculates = V x Ay{rk) by using Eq. ([2]). The sign in 
Eq. (j46p is taken for electrons with charge e = — |e|. It is chosen to minimize the 
energy [32j. For holes with charge e = |e|, we need to change the sign for v in Eqs. 
([2|) and ([TT]) . then Eq. (|36|) and the expectation value for energy (see below Eq. 
(jSOp ) retain the sign. 

If we take Sz = ^/2 and take into account that length unit is L, so 5^"^^ (r) should 
be replaced by 5^'^\r)/L'^, then 

^ N *.2 

^E^--4 = — E ^''Hr,-r,). (48) 
fc=i '^'^ Kj{k^j) 

The calculation of the expectation value, Eq. ([3]), with the Hamiltonian, Eq. 
(j48p . and wave function ^'T(ii), Eq. ([6]), gives (in n^/{ML'^) units) 

-7TU J ^t{R) <5(') (r-fe - r-) vl,^ (i?) di? = " 1) , (49) 

The total expectation value for the ground state energy including all terms of 
the Hamiltonian is 

E = ^ + ^ + NM L L , (50) 
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where the quantity Ai contains the Coulomb interaction and jelhum background 
term (see previous section). 

As in previous sections we have used the expression = 1 + i'{N — 1) for 
J\f, which is result of the cut-off parameter regularization [T5] of the logarithmic 
divergence, revealing in the bosonic representation, when two particles come up to 
each other. Substituting this expression for TV into expectation value for energy 
E, Eq. (j50p . we see the exact cancellation of terms containing the v factor. This 
can be also achieved putting in the formula for ground state energy = 0, i.e. the 
case of bosons. As the energy of bosons is lower than one for fermions or anyons 
with ^ 0, there appears a coupling of spin with statistical magnetic field for every 
particle or bosonization of 2D fermions and anyons. Of course, this effect occurs 
when the particles have the spin h/2. 

As next step of this coupling scenario might be the fluctuations of spins of bosons 
around of statistical magnetic field. Hence, the bosons might reveal an effective 
spin and look like as Fermi particles. However, Fermi gases with different spins are 
independent [33]. Thus, we have presumed that in the zero order approximation 
the spins of bosons interact only with each other and do not interact with spins 
of another fermions if they exist in the system. Interaction of spins might take 
place inside of some length- spin correlation radius, which for temperature T = 
we denote as ^o- The final feature of this model might be a destruction of spin and 
magnetic field coupling at the increasing of spin fluctuations when bosons become 
the fermions. This would occur when the gain in the energy of bosons due to 
fluctuations is equal to energy difference between the Fermi and Bose ground states. 

The interaction term of boson spins in the Hamiltonian we brought in the form 

N 

e-ro/io ^ g^^^ . 4 . (51) 
fc=l 

Here tq is the mean distance between particles. We introduced in this expression 
a factor e~^°/^°, which takes into account the exchange character of interaction of 
spins [34j- However, if the typical scale of .^o is the nearest interatomic distance, 
for our sreened by magnetic field spins, is to be assumed phenomenological and 
taken from experiment. 

For the T = following argument allowed us to establish the explicit 

form of Eq. (jSlh . The boson ground state energy was obtained when the term, 
Eq. (j46p . was included in the Hamiltonian. To get again the fermion ground state 
energy, we need to cancel it. Therefore, for dense (ro < ^o) boson gas, there should 
be Sk+s = —fibk/M. Comparing with Eq. (j47p we see that the summation index 6 
includes all — 1 numbers of index j. 

The expression for the Hamiltonian of bosonized infinite Coulomb anyon gas 
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with interaction of spins is 



1 

~ 2M 



N r 



\2 1 1 

Pk + A^in)) +MV{fk) +7;Y.T^ + 



N 



M 




k=l '- 



k=l 



(52) 



We can bring here the results of Ref. |16j and previous section. The ground state 
energy of the Coulomb anyon gas at > 2 with the high accuracy (see Figure 2) is 
described by formula, Eq. (|44p . with the expression for function /{v, r^), Eq. ()45p . 
For the HTCS, actual for us is the region > 9 (see below Figure 4 ), where one 
can put f{iy, rs) ~ f{i' = 0, r^). The ground state energy for the 2D Bose and Fermi 
gases is obtained putting in Eq. ()44p = and v = 1, respectively. 

The analogical to Ref. [T6] and Section 3 calculation, with the Hamiltonian, 
Eq. (f52l) . gave for > 2 the same expression, Eq. (l44|l . for ground state energy of 
bosonized anyons, but in it one needs to replace v by i/e~'"^/'»" (now is expressed 
in Ob units ). Considering below the bosonization of 2D fermions we have putted 
now V = \. 

To become the pure fermions bosons had to overcome the energy difference 



One needs a special remark here. Our previous calculations have been related 
to ground state of spinless or fully spin polarized fermions (electrons). It would be 
preferable to deal with one of normal, i.e. with no spin polarized, electron liquid 
state. However, the difference of their ground state energy is essentially lower |31] 
than accuracy of our calculations. 

We have applied the model for possible clarification of summarized experimental 
phase diagram of hole doped High-Tc superconductors, which was recently proposed 
by Tallon and Loram [U]. As it was said above it is believed that the main con- 
tribution into HTCS is provided by ab planes of set of Cu02 atoms and therefore, 
2D is responsible for it. The SC state occurs when holes transit into BEC. We note 
that for the ideal 2D systems BEC can exist only at T = [35j. At T / the 
evolution of fluctuations for the order parameter destroys BEC. However, if 2D gas 
is situated inside of third dimension this evolution is suppressed (with respect to SC 
see book of Abrikosov t36j). Furthermore, we used in the Hamiltonian, Eq. (|52p . a 
three dimensional Coulomb potential. Therefore, for HTCS the 2D BEC of bosons 
can also exist at not zero temperature. 

We used as gap for SC and studied its asymptoticses as function of density 
of holes Us = l/{'Kr1). Below we show that ~ p, where p, called in the experiment 
as concentration of holes, is the fractional part of hole per atom Cu. At big values of 
Tg or small Ug one can neglect the exponential factor in Eq. (j53p and A^ ^ Ug ^ p. 
At small or big n^, without this factor, A^ would have A^ ~ n^J'^ ^ p^/^ , but 
Q-'Ts/io suppresses this dependence to zero and the law of it depends from function 



7(1 -e-'-^/«o)/4/3(0,r,) 



(53) 
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^o(p)- At this limit of we presumed that the function Aq (p) coincides with 
experimental dependence Eg{p). Extrapolating this asymptotics of A^(p) to small 
values of p and equating it to Eg (p) one finds the empirical dependence (p) • For 

it Up) ~ i/Eg{p). 

The experimental values for energy are expressed in the Kelvin temperature, K, 
units. To be sure in the correctness of our 2D density of holes, we expressed it in 
cm~'^ units and compared with experiment. For elementary structural cell of almost 
all cuprate ab planes a = 3.81A and b = 3.89A. Assuming that it has one atom of 
Cu with p fraction of hole, the density is Uab = ^ab " P cm~'^, where Nab is number 
of elementary cells per Icm^ square. We have n^fe = 6.7472 • p ■ 10^^ cm~'^, from 
where p. At the optimal doping of holes, p ~ 0.16, we compared nab with 

experimental one of |37j . In it for Y-123 compound the square of three dimensional 
plasma frequency, being expressed in square of wave number units, for SC carriers 
has a value 1.1 • 10* cm~^. We find the density of holes per 1 cm^ and then, 
assuming that for Y-123 compound there are two Cu02 planes [5] per elementary 
cell, responsible for SC, and in the c axis the spacing is c = 11.7A, transform it 
into 2D n'^^ = 0.913710^^ cm~'^. Our optimal doping value for nab is 1.0795 • 
10^^ cm~'^. Therefore, typical density for SC carriers for the optimal doping state 
is ~ 10^^ cm^"^. From the experiment [37] also leads important information about 
approximate equality of SC carriers density to whole one. Thus, we can bring nab 
as SC density. A comparison of this value for nab with npQUE ~ 10^^ -i- lO^'^ cm~'^ 
for Fractional Quantum Hall Effect (FQHE) experiment [38j shows importance of 
the pure, not screened, Coulomb potential interaction, which is essence of FQHE 
[28], for HTCS. The screened Coulomb potential is supposed to be as justification 
[39j for the treatments based on the Hubbard model. 

On the Figure 3 we presented the summarized curve for pseudogap boundary 
energy Eg (Fig. 11 from paper [H]), as result of interpolation of sets of experimental 
data, SC gap energy Aq = 4:Kb Tc, which was evaluated by empirical formula 
Tc = Tc,max[l - 82.6(p - 0.16)2] with Tc^max = 95 K for Bi 2212 compound [40J, 
and calculated by formula Eq. (1530 SC gap energy for bosons A^ as function of 
concentration of holes p. Despite a A^(p) dependence has the qualitatively different 
form than experimental one, it is in accordance with the conclusion of Tallon and 
Loram that Eg, being going to zero at the critical doping concentrate pc ~ 0.19, 
separates BEC into two parts, where in first one the density of SC providing pairs is 
small (even though that conclusion is not expressed explicitly by Tallon and Loram, 
however, Fig. 10 of their paper more than obviously displays it). As we see from 
the figure the magnitudes of this and experimental diagrams have the same order. 
We would like to note that, apparently, it is the first estimate for BEC energy of 
canonical bosons, which can be appropriate for experiments of SC. 

Figure 4 demonstrates the p dependencies of and ^o- The spin correlation 
radius becomes to be sharply increased when p approaches pc, which might mean 
the vicinity of phase transition, where many and many spins of bosons are involved 
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in correlation process before the transforming of them into fermions. We estab- 
lished a correspondence with second important conclusion of Tallon and Loram that 
"the pseudogap is intimately connected with (though not equivalent to) short-range 
Anti Ferromagnetic (AF) correlations, which disappear at the same critical doping 
state" (pc), and results shown in Figure 4. The paper [13] shows that experimental 
short-range AF correlations scale like Eg{p) dependence and vanish at pc (see Fig. 
6 of this paper). It was also suggested in Ref. [H] that two these quantities are 
the temperature independent. Our treatment is not for the temperature T ^ 0. 
However, last remark allow us to do the connection between short-range AF corre- 
lations and correlation radius ^o- The presumption about correlations of spins of 
bosons would mean that there would exist the competing of these correlations with 
AF correlations, where nature of last ones is the spin interaction of closely situated 
fermions. Increasing of correlation radius would lead to revealing the Fermi like 
spin correlations of bosons in the extended region of sample and thus, assuming 
that in the first order approximation these spins now interact with ones of fermions, 
to suppressing of short-range AF correlations inside of this region. From this a 
disappearing of AF correlations at pc, where goes to infinity, would be natural 
and might be considered as indirect indication of destruction of BEG from canoni- 
cal bosons. Mathematically the short-range AF correlations will be proportional to 
l/(,o{p), because Eg{p) ~ l/(,o{p)- One might suppose that the interaction of spins 
of fermions with statistical magnetic field (as well as the statistical magnetic field 
itself) of anyons would have a hidden character and might not be revealed obviously 
in the experiment. In this case, the described experimental behaviour of short-range 
AF correlations might be also the implicit indication of existing of the statistical 
magnetic field. 

Next conclusion of Tallon and Loram paper [T3] is an independence of pseudo- 
gap boundary energy Eg and phase diagram for BEC. As we said above, in our 
treatment, without Eg, Aq would have a proportionality ~ p^/^. The phe- 
nomenological including of Eg determines the maximum Tc value for BEC phase 
and critical concentration pc- Close to pc Eg{p) merges with A^(p) determining the 
asymptotics for latter one. In our model Eg characterizes the cuprate material. As 
for almost all High- Tc superconducting samples the numerical values for cell spacing 
constants a,b are the same [5], according to our approach we should been obtain 
the similar Tc^max for all of them. The value for Tc^max depends also from number 
of Cu02 planes in the elementary cell |5j. We have considered the sample with one 
plane in the hoping that Eg{p) takes effectively into account this factor. 

Here we would like to say some words about the possible scenario for the phase 
diagram of HTCS, which might be derived from this model. Due to vicinity of the 
structural phase transition to superconducting state in the cuprate materials at the 
concentration of holes below optimal [5], we might do a hypothesis that the induced 
by it the mechanical strain would strengthen a quadratic striction and therefore, 
the phase transition into BEC would be not a second order, as it should be, but 
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first, close to second one [41j- For this case a pseudogap regime would find natural 
interpretation. It would correspond to meta stable phase of bosons. At pseudogap 
boundary, Eg, bosons would finally undergo a phase transition into fermions. The 
critical concentration, pc, at which the pseudogap energy. Eg, is zero, would be 
considered as critical point of first order phase transition. The first order phase 
transition scenario is possibly consistent with recent experimental observations of 
coexistence of pseudogap and superconducting phases described in the review of 
Pines [l2], because typically the coexistence of equilibrium phases exists in this 
order of phase transitions [l3]. In our case, these phases would be the small islands 
of meta stable phase of bosons and BEC situated in close area to BEC boundary 
in the phase diagram. For electron doped materials one may suppose that there 
is no structural phase transition. Therefore, the phase transition into BEC would 
have a pure second order, which is possibly seen in the experiments. At last, the 
hypothesis might clarify the possible independence of pseduogap energy. Eg, and 
BEC suggested in the paper [Mj. As it was said above the boundary. Eg, would be 
connected with mechanical strain property of cuprate materials, while condensation 
into BEC with coupling of spins of fermions with statistical magnetic field. However, 
the spin correlations of bosons, which are to be governed by Eg and suppress the 
BEC, might implicitly be relate two these quantities. 

Now about the spectrum of collective excitations, Eq. ([35]) . This spectrum 
has no gap at small wave vectors k and therefore, can not be responsible for SC. 
SC is provided by particles in the BEC. On the other hand, one might expect that 
the bosons with effective spins inside of correlation radius ^o(p) a-nd obeying the 
relation, Eq. ()35p . at big wave vectors k are not the Fermi liquid quasiparticles 
like. 

We considered the effect of bosonization of 2D fermions for the instance of HTCS. 
One might presume that it would be important for any 2D fermion gas. FQHE, gas 
of isotopes He 3 and many other objects might be affected by this effect. 

5 Summary 

We have introduced in the Hamiltonian of anyon gas the Zeeman term of the inter- 
action of spins Sz = h/2 of particles with magnetic field induced by anyon vector 
potential, i.e. statistical magnetic field. A calculation of the expectation value for 
ground state energy in the framework of variational approach with cut-off parame- 
ter regularization has exhibited the cancellation of terms connected with fractional 
statistics, which might be mean the bosonization of anyons due to a coupling of 
their spins with statistical magnetic field. For fermions, as the particular case of 
anyons, we have applied it for the possible clarification of summarized phase dia- 
gram of HTCS suggested by Tallon and Loram for the region below preudogap Eg. 
Additionally to Zeeman term we have phenomenologically introduced in Hamilto- 
nian the term, which was responsible for the correlations of effective spins of bosons. 
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Figure 3: The experimental pseudogap energy Eg, SC gap energy Aq = 4:Kb Tc 
(experiment for Bi 2212 compound), and calculated by formula Eq. ()53p SC gap 
energy for bosons in Kelvin temperature (K) units as function of concentration 
of holes p. 
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Figure 4: The mean distance between holes and spin correlation radius 
Bohr radius as units as function of concentration of holes p. 
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and connected them with Eg. The obtained phase diagram has been quantitatively 
close to experimental one, while model qualitatively described conclusions made in 
paper of Tallon and Loram. 
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